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Abstract 

In this study, a phenomenological model of the radio frequency (RF) behavior of a 
superconducting cavity fundamental power coupler is proposed by analyzing the simulation results 
of a transient beam loading process in an extremely over-coupled superconducting cavity. Using 
this phenomenological model, the calculation of the transient reflected power from a 
superconducting cavity under beam loading can be mathematically simplified to algebraic 
operations without solving the differential equation governing the transient beam loading process, 
while maintaining the calculation accuracy. Moreover, this phenomenological model can facilitate 
an intuitive understanding of the significant surge in the time evolution of reflected power from a 
superconducting cavity in certain beam-loading processes. The validity of this phenomenological 
model was carefully examined in various beam loading processes and cavity conditions, and the 
method based on this phenomenological model was utilized in the transient RF analysis of the 
superconducting cavity system of the CAFe Linac, achieving satisfactory results. 
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1 Introduction 


The calculation of the transient power reflected from a high-Q-value resonator under extreme 
over-coupling has practical significance in the field of microwave engineering because the excessive 
reflected power from this resonator under certain conditions may damage the radio frequency (RF) 
system. A superconducting RF cavity l! , adopted to accelerate the charged particle beams in an 


46] is a typical high-O-value resonator with intrinsic Q-values typically > 10° and 


accelerator ! 
loaded Q-values l usually > 10° due to the extremely low RF-power loss in the cavity wall owing 
to the superconductivity of niobium | °! at ultra-low temperatures |'°!, Unlike a normal conducting 
cavity with a loaded Q-value, typically in the order of 10° and near critical coupling, the transient 
power reflected from a superconducting cavity in certain beam-loading processes can attain values 
several times larger than the RF generator forward power l!!! , thus has a high potential to damage 
the RF power system, particularly the circulator, which is usually designed to withstand as large as 
the forward generator power considering a full reflection situation. Therefore, it is critical to 
calculate the time evolution of the reflected power from a superconducting cavity during various 
transient beam loading processes to provide the necessary information for the design of an RF power 
system and a machine protection system !'?!, the latter one can prevent the RF power system of an 


accelerator from being accidentally damaged and increase its availability. In existing calculation 
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methods |'?: 141, the transient power reflected from a superconducting cavity under beam loading can 
only be obtained by analytically or numerically solving the differential equation governing the 
transient beam-loading process inside the cavity, which is relatively complicated and can be time- 
consuming. For certain specific requirements, such as designing a high-level RF system concerning 
the RF power waveguide or the circulator between an RF power generator and a superconducting 
cavity, the transient power reflected from the superconducting cavity in various processes is the RF 
quantity of interest, instead of the time evolution of the RF fields inside the cavity !!>-!7], the analysis 


18-20] Therefore, a 


for the latter one is usually performed by low-level RF system engineering ! 
method to directly calculate the reflected power without numerically solving the beam-loading 
differential equations was researched. In this study, we reported a method based on a novel 
phenomenological model of a fundamental power coupler to fulfill the aforementioned goal, which 
has been facilitating the transient RF analysis of the superconducting cavity system of the CAFe 
Linac ?!3]. Thus, a phenomenological model of RF behavior of the fundamental power coupler 74! 
for a superconducting cavity under heavy beam loading was obtained by investigating the simulation 
results of a superconducting cavity in a transient beam-loading process without directly analyzing 
the corresponding differential equations governing such a process. Heavy beam loading refers to 
the case where the beam power is several orders of magnitude larger than the RF power loss in the 
superconducting cavity wall. Using this phenomenological model, the transient power reflected 
from a superconducting cavity can be converted into algebraic operations without solving the related 
differential equations. To verify its validity, the transient reflected power in several beam-loading 
processes was evaluated using the proposed model, and the results are compared with those obtained 
with the numerical method needing to solve the related differential equations. Moreover, the 
phenomenological model can analytically deduce the maximum transient reflected power during a 
superconducting cavity failure recovery process in an extremely high-availability accelerator, such 
as the one used for the accelerator-driven sub-critical system (ADS) application. 


2 Existing calculation method and result analysis 
2.1 Mathematical model on transient beam loading 

In the conventional calculation method, it is necessary to solve the following equations 
governing the transient beam-loading process analytically or numerically to obtain the transient 


power reflected P, (t ) from a superconducting cavity under the Gaussian bunched beam loading 


[25]; 
oy iAw)V_ osř Ft) io - fs 
ETs a taoT O= i fa PO goan (i)Tye } (1) 
and 
P(i)=P (=r ()r Trao] dV, (t) (2) 


where i is the imaginary unit, V.(t) =V, (t)exp[i- (0) | is the phasor of the cavity voltage 


without the transit time factor Tp; @, @,;,and A@=@-@, are the angular frequency, loaded 
half bandwidth, and detuning of the cavity, respectively. ~, and f / Q are the intrinsic shunt 
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impedance and r over Q of the cavity without Ts, respectively; is the coupling factor, P, (t) 
is the forward generator power, O(t) is the phase of the effective driven current from the generator, 
I,)(t) is the time-averaged DC beam current, and ¢,;=@0o, is the half phase width of the 
Gaussian bunch (the effect of the beam current intensity on the bunch length is neglected here). 

It can be mathematically proven that V(t) can be analytically solved using the Laplace 
transformation method when the time evolutions of P, (t) and 7 w(t) in Eq. (1) conform to Eqs. 
(3) and (4) respectively, which is sufficient to describe a typical beam-loading process in a 


superconducting cavity >); 


Pı 0<t<t 
3(t-4,) 2(t-4,) 
P+ E T (P Pa) t <t<t,+At, 
1 1 
P (t)= P, t, +A <t<t, , (3) 
3(t-t,) 2(t-4) 
P, l T 2 T = |(P.,-P.,) t, <t<t,+At, 
2 2 
Ps t2t,+AT 
Ly, O<t<t,, 
and I,,(t)= 1 Ty ba Stha (4) 
Iy to z t< ty3 


Substituting the resultant V, (t) into Eq. (2), P(t) can therefore be analytically obtained . 
2.2 Simulation result and analysis 


The designed RF parameters of the HWR010!°! and RF power generator of CAFe Linac, along 
with the designed beam parameters, are listed in Table 1, the transient beam loading process can be 


simulated using these parameters and Eqs. (1)-(4). 


Table 1. RF and Beam Parameters for the evaluation 


Parameters Value 
Mode frequency, fo (MHz) 162.5 
r over Q (Q) 147.84 
Intrinsic quality factor, Qo 8x108 
Optimum coupling, Bopt 800 
Transit time factor, Tir 0.8 
Forward generator power, Pg (kW) 6.7 
Phase of the generator current, 0 -7/6 
Synchronous phase, ¢ -7/6 
Half phase width of the bunch, Ado.s m/12 
Design beam current, /b (mA) 10 


Generator power rising edge, Afi (us) 10 


Generator power falling edge, Ati (us) 10 


The time evolution of the generator power and beam current in this process are shown in Fig. 
1 (a). At f =0, the generator is powered up at 6.7 kW, and the cavity field buildup process 
begins. At f, =AT , the continuous-wave (CW) proton beam began injecting at the designed 
current of 10 mA as the cavity field reached a certain steady state I with an amplitude of 1.72 MV 
and phase of zero. Owing to the beam loading effects, the field in the cavity evolved to steady 
state II with an amplitude of 1.00 MV and phase of —29.8°. At t, =2AT , the generator power 
was suddenly shut down and the beam continued passing through the cavity at 10 mA current, 
which is typical for a Linac used for ADS application to avoid a thermal shock to the reactor core 
due to the beam trip. This resulted in steady state III with a cavity field amplitude of 0.987 MV 
and phase of -150°. At t, =3AT , the generator power was recovered to 6.7 kW, and the field in 
the cavity restored to steady state II. During the entire process, the detuning angle of the cavity Y 
maintained its optimum value for stage II with the designed beam current intensity, which is as 


follows: 
Y opt Ñ —ġ ~ 1/6 T (5) 


In the simulation, we used AT x8r to ensure that the cavity fields can reach the 
corresponding steady state at each stage. Fig. 1 (b)-(f) presents the time evolution of the trace of 
the cavity voltage on the complex plane, cavity voltage amplitude, phase, reflected power, and 
kinetic energy gain of the central particle in a beam bunch, respectively. The simulation results 
shown in Fig. 1 were analyzed to obtain a phenomenological model of the fundamental power 


coupler for a superconducting cavity. 


202306.00050v2 


chinaXiv 


—P (kF) 
g 
F- |— (ma) 


0 0.01 0.02 0.03 0.04 0.05 0.06 Be E , 15 l 23 
(MV) 
t(s) r 
(a) (b) 


S 
č 
0 0.01 0.02 0.03 0.04 0.05 0.06 0 0.01 0.02 0.03 0.04 0.05 0.06 
t(s) t (s) 
(c) (d) 
1.24 
1] 
0.84 
œo 064 
S z ] 
A = 2] 
ie y“ 0.24 
o4 
-0.24 
-0.44 
-0.64 
t t t + t - r - - 7 t v t d 
0 0.01 0.02 0.03 0.04 0.05 0.06 0 0.01 0.02 0.03 0.04 0.05 0.06 
t (s) t (s) 
(e) (f) 


Fig. 1 (Color online) A superconducting cavity transient beam-loading process 
(a) Generator Power and Beam Current (b) Trace of the Cavity Voltage (c) Amplitude of the Cavity Voltage 
(d) Phase of the Cavity Voltage (e) Reflected Power of the Cavity (f) Energy Gain of the Central Particle 
Fig. 1 (c) demonstrates that the amplitude of the cavity voltage in steady state III is the same 
as that at steady state II. 


Va| = oa © 
At steady state II, the cavity voltage Fe is the superposition of the generator-induced 
voltage va and the beam-induced voltage Va : 
Vo = Var o (7) 


In steady-state III, owing to the absence of the generator power, the cavity voltage V; is 


solely induced by the beam: 


Fa =h (8) 


c 


Using Eqs. (6) — (8), the following is obtained: 


Kal 


T+T]. ©) 


Note, the beam current intensity remains constant during stages II and MI, V,a = V5 =F, 


and Eq. (9) can be rewritten as follows: 


KF 


AE (10) 


This indicates that the amplitude of the beam-inducing voltage is the same as that of the 
cavity voltage produced by the generator along with the beam, which is the case for a 
superconducting cavity under optimum coupling and detuning |"). 

As shown in Fig. 1 (d), although the cavity voltage amplitudes at steady states II and III are 
identical, there is a difference of 27 / 3 in their phases, which results in the deceleration of the 
beam at steady state III, as shown in Fig. 1 (f). This can be explained through the conservation of 
energy principle, where the field in the cavity is induced solely by the beam in this state. In 
addition, Fig. 1 (e) indicates that the reflected powers at steady states II and III are significantly 
different. At steady state II, the optimum detuning and coupling of the cavity ensure zero-power 
reflection, whereas at steady state III, the cavity functions as a klystron, where the beam power is 
drawn out by the cavity and radiates to the waveguide through the coupler. Note, the transient 
reflected power is three times as large as the generator forward power at t, =3AT when the 
generator power is recovered. The excessive reflection power is due to the change in the stored 


energy inside the cavity, which corresponds to the last term in Eq. (2). 


3 Relationship between the voltage waves in the waveguide 
3.1 Amplitude relationship between the voltage waves 

According to the terminology of microwave engineering P7], the reflected power from the 
cavity P(t) in the previous sections is termed as “backward power in the waveguide” and is 


referred to as “backward power” for conciseness in the following sections. Meanwhile, to maintain 


brevity, the symbol a is used to denote the amplitude of a phasor a, thatis a= la 


According to the theory of transmission lines, the traveling wave power in a waveguide can be 
expressed by the amplitude of the normalized equivalent voltage wave v as follows 78): 
Pav , (11) 
where v= v/ JZ, , Vis the amplitude of the normalized equivalent voltage wave, and Z, is the 
characteristic impedance of the waveguide. 


Based on Eq. (11), the following can be obtained: 


v=4P . (12) 


In the following analysis, the normalized equivalent voltage wave amplitude v that was 
determined using Eq.(12) is termed “voltage wave” for brevity. 
The backward voltage wave (“backward wave”) v, in the waveguide comprises the 


following two components: 


1) Reflected voltage wave at the coupler (“reflected wave”) v, , which is the direct reflection 
of the forward voltage wave (“forward wave”) from the generator v, at the coupler. 
2) Radiated voltage wave from the cavity (“radiated wave”) v, , which corresponds to the 
power radiated out of the cavity back into the waveguide through the coupler. 
Ve HV +V, - (13) 


As shown in Fig. 1(e), when the cavity is under optimum coupling and optimum detuning, and 


the beam current is at the designed value, the backward wave at steady state II vanishes, i.e., 


os = 0; therefore, the following can be obtained: 


e,2 = V2 2 (14) 


This indicates that the reflected and radiated waves had the exact same amplitudes and opposite 
phases at steady state II. 
The amplitude of the cavity voltage at steady state III is equal to that at steady state II, note 


that v, €V, ; therefore the following can be obtained: 


e3 =V * (15) 


For a cavity with beam loading, the cavity voltage V, is the superposition of the generator- 


c 


induced voltage V, and the beam-induced voltage V, , i.e., V,=V,+V,. Similarly, the radiated 
wave v, through the coupler can also be composed of two parts, including the contribution from 


the generator Va and that from the beam v, : 


ee (16) 
with the following relationships: 
"he (17) 
v £V, 


At steady state III, the generator is powered down, and Vig = 0. Note that V3 a Veg ; therefore, 


V,a =0 and the resulting backward wave can be written as follows: 


SVa > (18) 


Fig. 1 (e) demonstrates that the backward power at steady state III is equal to that at steady 


state I, indicating that the amplitudes of the corresponding backward waves should be identical: 


Vok,3 = Vo + (19) 


Note, there is no beam in the cavity at steady state I due to the extreme over-coupling of the 
cavity (8 >>1), and the reflection coefficient of the cavity is I —1 ; therefore, the following can 


be obtained: 


a ae (20) 
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The forward power from the generator remains constant at steady states I and II: 


Veo = Ver: (21) 


Using Eqs. (14)-(21), the following amplitude relationship can be deduced: 


Viks T Yor TY, n2 T Vh T Vea TVer o (22) 


and thus 


Vio = Veo: (23) 


This indicates that the amplitude of the wave directly reflected from the forward wave at the 
coupler is exactly equal to the amplitude of the forward wave at steady state II. In simple terms, the 
forward wave is fully reflected by the coupler. 

The proportion of the forward wave reflected directly from the coupler is determined by the 
geometric structure of the coupler alone, which corresponds to the coupling factor B. Therefore, if 
P is fixed, the proportion of the forward wave directly reflected by a certain coupler should also be 
fixed. Note that B remains constant throughout the entire process; therefore, the following is 


obtained: 


Vai = Ve (i = 1..3) ; (24) 


This indicates that the forward wave is fully reflected at the coupler during the entire process. 
3.2 Phase relationship between the voltage waves 


Let the phase of Vio be 7 , and using Eq. (23), Vig can be expressed as follows: 


Vig = Ve + (25) 


Substituting it into Eq. (14) results in the following: 


Vor = Veo > (26) 


At steady state I, the cavity voltage V,, is induced by the generator alone. Considering the 


phase of the equivalent current of the beam as the reference phase, the phase of Vo is as follows: 
arg(V..)=O+y A (27) 
where arg(V,;) indicates obtaining the principal argument of V. : 


The phase of the cavity voltage V, at steady state II was the synchronous phase: 


arg(V..)=¢. (28) 


For a coupler with a fixed geometric structure, the phase advance of the electromagnetic wave 
transmitted through the coupler must also be fixed; thus, the phase difference between Ves and 
Ve 


el 


is equal to that between V,, and Fig i 


arg(v,,)-arg(v,;)=arg(V,,)—are(V,,)=¢-0-y ; (29) 


When the cavity is under optimum coupling and detuning, and the beam current is at the design 
value, the equivalent impedance of the beam-cavity system is real, indicating that the cavity voltage 


and the equivalent driving current from the generator should be in phase as follows: 
g=0. (30) 


Substituting this into Eq. (29) to obtain arg(v.z )-arg(,; ) =-y , and using Eq. (26), Ver 


can be expressed as follows: 


ev (31) 


Similarly, the phase difference at steady state II between Via and Veo is equal to that 


between V,, and Vee 
arg(v,.)-arg(v,;) = arg(V,.)-arg(V,.) : (32) 


At steady state II, the phases of the generator-induced voltage Fa and beam-induced voltage 
Vs are arg(V,2) =0+y and arg (V2) =a+y , respectively. Substituting them into Eq. (32) 


results in the following: 
arg(v,,)-arg(v,))=7-8 : (33) 


At steady state I, the radiated wave Va is entirely depends on the generator because the cavity 


voltage V,, is induced by the generator alone: 


Ver = Ver (34) 


Note that the forward power from the generator remains constant during stages I and II; 


therefore, the following can be obtained: 


a (35) 


g1 


The following can be obtained by using the two expressions above and Eq. (31): 


V> =v e” : (36) 


B 


The expressions above, Eq. (22), and Eq. (33) can be used to obtain the following: 


Va = vpe T (37) 


Using Eq. (12), Veo can be rewritten as follows: 


Wee eae (38) 


g,2 


Using the above expression with Eqs. (25), (26), and (35), the phasor relationship diagram 


of the various voltage waves in the waveguide can be obtained as shown in Fig. 2. 


Fig. 2 Phasor relationship diagram of the voltage waves in the waveguide 
Substituting Eqs. (26), (35) and (37) into Eq. (16) with the phase relationships at steady 
state II in Eqs. (5) and (30) to solve the resultant equation obtains v,, =2v,,cosy , which is then 


substituted into Eq. (36) to obtain the following equation: 


v., =2v,,coswe” . 39 
g.2 #2 


Using Eq. (12), the expression above can be rewritten as follows: 


v, =2,/P.,cosye” . (40) 
82 82 


4 Modeling the RF behavior of the fundamental power coupler 
4.1 Reflected wave and generator-induced radiated wave 
Generalization of Eqs. (25) and (40), assuming that they are always true for a fundamental 


power coupler in the case of extreme over-coupling, can be expressed as follows: 


vV, =-V; y (41) 
and Ve = 2,[P,cos ye” . (42) 


Eq. (41) indicates that the forward wave from the generator is always completely reflected by 


the coupler. Using Eq. (12), Eq. (41) can be rewritten as follows: 


es E (43) 


4.2 Beam-induced radiated wave 

To calculate the transient backward power in the waveguide, it is necessary to generalize 
Eq.(38). To simplify the problem, the case where the field in the cavity is induced solely by the 
beam was analyzed to obtain the general expression for the beam-induced radiated wave, because 
the field in the cavity satisfies the principle of linear superposition. 

When the field in the cavity is induced solely by the beam, the power radiating from the cavity 
through the coupler can be obtained using Eqs. (16) and (11) as follows: 


R= v = v? : (44) 
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Meanwhile, the corresponding coupling factor 4 can be expressed as follows: 


2 2 
pare ale ote (45) 
Pa K W 
where P, =V; /r. is the power loss in the cavity wall owing to the beam-induced field in the cavity. 


With Eq. (45), the following is obtained: 


v, © BV, . (46) 


According to the definition of the detuning angle, the beam-induced voltage V, should conform 
to the following: 
V, £ cosy . (47) 


To deduce the relationship between V, and £ , a transient beam-loading process is simulated 
with Eqs. (1) - (4) for different values of . At f =0, the CW beam injection begins and the 
field in the cavity is solely induced by the beam. The field in the cavity reaches a steady state after 
a sufficiently prolonged time and the detuning angle of the cavity remains constant throughout the 
entire process. The simulation result of the time evolution of the amplitude of the beam-induced 
voltage in this process is shown in Fig. 3, where Vs is the amplitude of the steady state beam-induced 


voltage corresponding to the coupling factor Bo. 
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Fig. 3 (Color online) Time evolution of the amplitude of the beam-induced voltage 


Based on Fig. 3, the steady-state value of V, conforms to V, œ B~. Using Eqs. (46) and (47) 


the following is obtained: 


(48) 


Note that v,, is for the optimum detuning and coupling; therefore, the general expression for 


V, is as follows: 


_ cosy Boos 


= ; 49 
Vy COSY 5, B b,2 (49) 


By assuming that the phase shift of the wave transmitted through the coupler remains constant 
irrespective of the coupling factor B, i.e., v, || as , and noting that cosy and cosy,,, are both 


real, B and Bop: are both real and positive, the following can be derived from the above relationship: 


= Cosy Pps — 


= . 50 
i cos Y opt p pe i ) 


Substituting Eq. (38) into the above expression results in the following: 


7 _ cosy Boy Peo alow). (51) 
COSY opt B 


Using the following relationships in an extreme over-coupling superconducting cavity !"41, 


Y opt Ñ —ġ 
h orh 
Bopi = P = y? 


ag . (52) 
Po ~F, =LoTre * V.cosd 


ae E 
B ii, (5o 


Eq. (51) can be converted to the following: 
— Ate r , 0 
Yaa * Ilo JO cospe) (53) 


4.3 Phenomenological model of the fundamental power coupler 
Using Eqs. (42), (43), and (53), the steady-state phenomenological model of the fundamental 
power coupler for a superconducting cavity with the Gaussian bunched beam loading are formulated 


as follows: 


Ve = 2,[P,cos ye” f (54) 
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Ah, 


v, 7 Tle è p2 cos ye) 


By replacing the beamform factor exp(—A¢,. / 2) in Eq. (53) with the corresponding 
equation, the phenomenological model of the coupler for other statistical distributions of the beam 
bunch can be readily obtained. For cases where the beam consists of short bunches, the beam form 
factor can be simply approximated as 1. In a case without beam loading, substituting /,,=0 into 
Eq. (53) obtains v, =0, and the phenomenological model of the coupler is further simplified. 

Using Eqs. (13), (16), and (11), the backward wave is obtained as follows: 
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Vig EV, EV +V, , (55) 


and the corresponding backward power is as follows: 


—* 


P, =y = Vix “Vik > (56) 


where Va "denotes the conjugate of a complex quantity vy . 

Note that Eqs. (3) and (4) requires P(t) and J,,(t) keeping constant for a short duration 
which is much larger than the cavity time constant after each time of their variations . 

IfP (t) and J,,(¢) can be approximated to be constant during a short duration after each time 
of their variations, and while accounting for the cavity detuning Aw and cavity filling process 
characterized by the cavity time constant T , a complex detuning attenuation factor A(t) can be 


obtained as follows: 
4 ken tl 
A()=exp| (i044) i (57) 
T 


and its counterpart, a complex detuning growth factor B(t) , can be derived as follows: 


B(e)=1-4(0) =1-ep (a04 t)r (58) 


The transient-state phenomenological model of the fundamental power coupler can be obtained 
by multiplying the corresponding factor to each expression in Eq. (54), with which the transient 
power reflected from the superconducting cavity in various beam-loading processes can be easily 
calculated using algebraic operations. 

Using the phenomenological model of the fundamental power coupler formulated by Eqs. (54) 
- (58), the calculation of the transient power reflected from a superconducting cavity is converted 
from the cavity-side of the generator-coupler-cavity system to the generator-side, thus eliminating 
the complexity involved in the time evolution of the cavity fields during the beam loading process. 


Therefore, this phenomenological model can significantly simplify the calculation. 


5 Verification and application of the phenomenological model 
5.1 Transient beam loading process A 


To validate the phenomenological model of the coupler, the time evolution of the reflected 
power from the superconducting cavity during the transient beam-loading process from Section 2 
was recalculated using the phenomenological model of the fundamental power coupler that was 
formulated using Eqs. (54) - (58). The following process exemplifies the procedure for calculating 
the transient reflected power from the cavity using this phenomenological model. 


Stage I: At ż=0 , the generator is powered up, and P, =P,,; since the beam is absent, 


gl? 


ILo =0. Substituting these values into Eq. (54), v,, and v,, can be obtain as follows: 


n=- . (59 
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and v,,(4)=0 . (60) 


The corresponding steady-state generator-induced radiated-wave Vega can also be obtained 
using Eq. (54). By multiplying it with the growth factor, Eq. (58), the transient generator-induced 


radiated wave Vet (t) can be obtained as follows: 


a ee ee di- exp -[iaeos agi ; (61) 


Stage II: At ż=AT , the forward power from the generator remains constant, that is, 


P, = P., = P, ; therefore, the following is obtained: 
val) =a =y > (62) 


and 


vya (t)= va (t) - (63) 

Beam injection begins with 7, ={,o2 ; by substituting this into Eq. (54), the corresponding 
steady-state beam-induced radiated wave Dine can be obtained. Multiplying this by the growth 
factor obtained by Eq. (58) with the time displacement of AT , the transient beam-induced radiated 


wave v(t) can be obtained as follows: 


Fal) Fas Berar) =f ex (aoste n|| (64) 


Stage III: At ¢=2AT, the generator is powered down, P, =0, and is substituted into Eq. 
(54) to obtain the following: 


mee (65) 


Multiplying Ven (t) in Eq. (63) using the attenuation factor in Eq. (57) with the time 
displacement of 2AT , the transient generator-induced radiated wave Ves (t) is obtained as 


follows: 


O=) Al 1-207) = (0)e09|-{io+4\(-207)] (66) 


Because the beam current intensity remains constant, the following can be obtained: 


va (t)= valt) 3 (67) 


Stage IV: At t=3AT , the generator is powered up again, and P, = P, . Substituting into Eq. 
(54) results in the following: 


< 


Veg Ba - (68) 


Using Eq. (54), the corresponding steady-state generator-induced radiated wave ea can 


also be obtained. Multiplying this by the growth factor, Eq. (58), with the time displacement of 
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3AT , the transient generator-induced radiated wave v,a (t) can be obtained (the contribution from 


v,3(¢) at this stage is negligible because it is sufficiently attenuated). 


rill) = Nga Blt 4-287) =e" ex (aost) an| ; (69) 
T 


Because the beam current intensity remains constant, the following is obtained: 


Td =l) - (70) 

The transient reflected power in each stage of P (t) (k=1.4) can be obtained by 
substituting the expression of each voltage wave in each stage of Eqs. (59) - Eq. (70) into Eq. (55) 
and further into Eq. (56). The time evolution of the reflected power from the cavity, P (t), during 
the entire process can be obtained by splicing P(t) according to the time sequence. By 
substituting the parameters listed in Table 1 into the resulting piecewise analytical expression for 
P (t), the resulting time evolution of P (t) is indicated by the yellow curve shown in Fig. 4. To 
verify this result, the time evolution of P (t) obtained using Eqs. (1}(4) is indicated by the red 
curves as a comparison. The thinner yellow curve is plotted above the thicker red curve; if the 
difference between the two results is sufficiently small, the yellow curve should be enclosed by the 
red curve, that is, the compound curve consisting of the red and yellow curves should have a closed 
red boundary. The red boundary of the compound curve is open only at one location where the sharp 


peak reaches its maximum; the relative difference between the two curves is only 1.91%. 
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Fig. 4 Transient reflected power in process A 

5.2 Transient beam loading process B 

To further verify the phenomenological model, the reflected power was calculated in a more 
complex transient beam-loading process using Eqs. (54)-(58) under the time evolution of the 
generator power and beam current intensity, which is shown in Fig. 5. The calculation verifies where 
the cavity largely detunes causing a complicated trace of the cavity voltage on the complex plane. 
However, under real-world conditions, the likelihood of such a transient beam-loading process 
occurring is a low. In contrast to the previous transient beam loading process A, the time evolution 


of the generator power and beam current in this process varied between two nonzero values. 
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Fig. 5 Time evolution of the generator power and beam current in process B 
The following calculation consists of two situations depending on the conditions of the cavity 
detuning and coupling. 


Cavity Condition I 
The cavity significantly deviates from its optimum detuning and coupling for the designed 


beam current in stage II, presenting a cavity detuning angle y = 2y,,, =7/3 and cavity coupling 
factor Z =42/32, ~ 653 during the entire beam loading process. At t = 0, the generator is 


powered up at 6.7 kW and the cavity field buildup process begins. At ¢,= AT , the CW proton 


beam at the designed current of 10 mA is injected as the field in the cavity reaches a certain steady 
state I with an amplitude of 1.10 MV and phase 30.6°. Owing to the beam loading effects, the field 
in the cavity evolves to the steady state II with an amplitude of 0.61 MV and phase of -5.28°. At 


t, =2AT , the generator power suddenly shuts down and the beam continues to pass through the 
cavity at a current of 3.33 mA, which is 1/3 of the designed beam current, resulting in steady state 
III of the cavity field with an amplitude of 0.417 MV and phase of 13.4°. At t,=3AT, the 


generator power is partially recovered to 5.36 kW, and the field in the cavity evolves to the final 
steady state IV with an amplitude of 0.790 MV and phase of 21.6°. The corresponding trace of the 
cavity voltage on the complex plane during the entire beam loading process is shown in Fig. 6 (a); 
the trace form is significantly different than that shown in Fig. 1 (b), where the cavity is under 
optimum detuning and coupling for the designed beam current at stage II. 
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Fig. 6 (Color online) Trace of the cavity voltage and transient reflected power: 


(a) and (b) demonstrate moderate detuning, whereas (c) and (d) demonstrate extreme detuning. 

For brevity, only the expression for Ves in phase III during this beam-loading process is given 
below as a demonstration. Stage III starts at ¢=2AT and v,,(¢), during which it can be viewed 
as a linear superposition of the following two sources: 

1) v,,(t) induced by P, and all the previous sources P,, for 1<2AT. 

2) v,;'(t) induced by the new source P, for t>2AT. 

At t=2AT, P , jumps to P, the sources of va (t) vanish, and v(t) in stage III is 


attenuated. 


va (t) =v, (2AT) -A(t >t-2AT). (71) 
Meanwhile, v,,'(¢) starts to rise with its driving sources P,,: 
vas (t)= Vos" B(t>t-2AT) , (72) 


where Ve ' corresponds to the steady-state value of Ves ' determined using Eq. (54). 


Using the above two expressions, Vis (t) can be obtained as follows: 


Vos (£) = Vaz (t)+ vas (0) oo 


=v,)(2AT)- A(t >t-2AT)+v o3" B(t>t-2AT) 


(73) 


The expression for Vya (t) can be obtained in a similar manner. 

By substituting the parameters listed in Table 1 to evaluate P (t), the resulting curve is shown 
in yellow in Fig. 6 (b), whereas the red curve is obtained by directly solving Eqs. (1) - (4). The 
figure demonstrates that the red boundary of the compound curve is closed during the entire beam- 


loading process, implying that the transient reflected powers from the cavity obtained by the two 
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methods show sufficient agreement. 


Cavity Condition I 

The cavity is at extreme detuning with the detuning angle w = 72/15, which approaches 1/2 
and significantly deviates from the optimum coupling for the designed beam current at stage IT with 
the cavity coupling factor 2 = f2/3Bops = 653 during the entire beam-loading process. Owing to 
the extreme detuning of the cavity, the corresponding trace of the cavity voltage on the complex 
plane in this process becomes significantly more complicated than that in all the previous situations, 
as shown in Fig. 6 (c). By substituting the parameters listed in Table 1 to evaluate P (t) , the resulting 
curve is shown in yellow in Fig. 6 (d), whereas the red curve is obtained via the direct solution of 
Eqs. (1) - (4). The red boundary of the compound curve is closed during the entire process, 
indicating that the transient reflected power from the cavity obtained by the two methods is in 
sufficient agreement, despite extreme cavity detuning. Note that the commonly used conventional 
method for numerically solving Eqs. (1) and (2) would be significantly time-consuming when 
dealing with the beam-loading process with a cavity detuning angle that is approaching 7/2, as in 
this case. 

Additional comparisons for various beam-loading processes and cavity conditions were 
conducted between the transient reflected power obtained from the phenomenological model using 
the formulated Eqs. (54)-(58) and that obtained by numerically solving Eqs. (1) and (2); the 
results were in sufficient conformance with each other as long as the cavity was extremely over- 


coupled. 
5.3 Maximum reflected power during the cavity failure recovery process 


As shown in Fig. 2, the peak of the transient reflected power occurs when the generator power 
is recovered while the beam continues to pass through the cavity during the period without the RF 
power feed from the generator, which is a realistic situation that occurs in an accelerator requiring 
an ultra-high beam availability, such as that in the CAFe Linac. The generator power driving the 
superconducting cavity is instantly cut off by the machine protection system to protect the cavity 
when the cavity quenches; thus, the trend of quenching is presented, or the instability of the cavity 
fields is detected. The generator power is restored after the problem of the cavity is resolved while 
the beam continues to pass through the cavity during the entire process. A technique for cavity 
failure recovery is currently under development for the CAFe Linac to ensure that its beam 
availability can eventually meet the ultrahigh demand of beam availability for the ADS application. 
Therefore, for such accelerators, it is critical to calculate the maximum reflected power during a 
cavity failure recovery process to provide indispensable information for the design of relevant parts 
of an RF system. By applying the phenomenological model formulated by Eqs. (54)-(58), this 
maximum value can be deduced without the specific parameters of the accelerator. 

At the moment when the generator power is recovered, V, = 0, and based on Eq. (17), Ve =0; 


therefore, the backward wave can be expressed using Eqs. (54) and (55) as follows: 
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Ag 
wr AR aR thate (0, cospe oa 
= -/P, + ier) 

and the amplitude of v is as follows: 

v= LJP +ye < Ria +|- (75) 


The equality expressed in Eq. (75) is only true if —/P, and ve") are in phase, where 
z-O+Ww=n (76) 


or 


y=0. (17) 


For a superconducting cavity accelerating CW or long macro—pulse beam, the cavity is usually 
optimally detuned to minimize the required RF power from the generator, and using the first 


equation in Eq.(52), the following is obtained: 
V= =O . (18) 
Eqs. (77) and (78) result in the following: 


Vor =0=0 , (79) 


which represents the on-crest acceleration of an ultra-relativistic beam. 


Furthermore, the maximum amplitude of wy is as follows: 


Vok,max = KA E A = JP +v, . (80) 


Using Eq. (22) to obtain v, =v, =,/P, under an optimum detuning and coupling with the 


designed beam current, Va max Can be further simplified to the following: 
Vok,max = 2 JP, 3 (8 1) 


and the corresponding maximum backward power, or the reflected power from the cavity, 


P. is written as follows: 


r,max 9 


P. = Veri Z 4P, $ (82) 


r,max 


This indicates that the upper limit of the transient reflected power during the cavity failure 
recovery process is four times that of the forward generator power for a superconducting cavity 
under optimum detuning and coupling to the acceleration of the beam at the designed current. 
Therefore, the design of the relevant parts of an RF system should incorporate a transient power 
surge five times (including a one-time contribution from the forward generator power) that of the 


forward generator power to ensure the safety of the instrument under such conditions. 


6 Conclusion 


Using the phenomenological model of the RF behavior of the fundamental power coupler 


19 


formulated using Eqs. (54)-(58), the mathematical calculations for the evaluation of the transient 
power reflected from a superconducting cavity in the presence of beam loading can eventually be 
converted to algebraic operations without solving the relevant differential equation numerically. 
This can be done while maintaining the calculation accuracy, which signifies the vast simplification 
of the conventional calculation process and significant reduction in the computational time, 
especially in the early design stage of parameter optimization. Moreover, for the beam-loading 
process, the difference between the theoretical results predicted by solving Eqs. (1)-(2) and the 
corresponding experimental values can typically reach 5%, which is larger than the difference in the 
results from the phenomenological model and directly solving the corresponding differential 
equation, indicating that the phenomenological model is sufficient to predict the transient reflected 
power in such a process. 

Meanwhile, the simple and vivid physical demonstration provided by this phenomenological 
model can facilitate an intuitive understanding of the time evolution of the reflected power from the 
superconducting cavity in various transient beam-loading processes, especially the significant surge 
in the reflected power during a certain beam-loading process, such as that simulated in Sect. 2. To a 
certain extent, the phenomenological model can be regarded as an extension of the classical 
transmission line theory of microwave engineering, where a radiated wave emerging from an 
extremely over-coupled active microwave impedance is introduced under the circumstance of the 
superconducting acceleration. 

Note that by using a simple beam-loading process to deduce the phenomenological model, the 
model may appear to be valid only for the beam-loading process with a simple time evolution of the 
generator power and beam-current intensity. However, the model can be validated for any transient 
beam-loading process because the time evolution of excitations described by an arbitrary continuous 
curve can always be approximately described by a staircase-like polygonal curve similar to Eq. (3), 
where a piecewise function with more nodes is required to approximate the original curve for the 
time evolution of a more complex excitation. For each stage of the piecewise function, the duration 
needn’t to be significantly longer than the cavity time constant to reach the corresponding steady 
state, instead, it can be significantly shorter than the cavity time constant to describe a notably more 
rapid process. 

The method reported in this study based on the phenomenological model of the RF behavior 
of the fundamental power coupler was successfully utilized in the transient RF analysis of the 
superconducting cavity system of the CAFe Linac over the past several years and has achieved 
satisfactory results, while vastly simplifying the calculation process and significantly reducing the 


computational time. 
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